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Abstract 

We discuss two simple variational approaches to quantum wells. The trial harmonic functions 
analyzed in an earlier paper give reasonable results for all well depths and are particularly suitable 
for deep wells. On the other hand, the exponential functions proposed here are preferable for 
shallow wells. We compare the shallow-well expansions for both kind of functions and show that 
they do not exhibit the cubic term appearing in the exact series. It is also shown that the deep-well 
expansion for the harmonic functions agree with the first terms of perturbation theory. 
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I. INTRODUCTION 



In a recent paper we discussed the application of the variational method to a Gaussian 
well.- We showed that harmonic variational functions yield reasonably accurate results for all 
well depths despite the fact that the shallow-well expansion for the ground-state variational 
energy does not agree with the exact expansion beyond the leading term. 

Students at introductory courses on quantum mechanics are commonly encouraged to 
solve the Schrodinger equation for a particle in a one- dimensional square box with finite 
walls. They thus learn that the exact solution decays exponentially as \x\ — > oo. After such 
an experience they would find it unreasonable the choice of harmonic variational functions 
that do not exhibit the behavior at infinity expected for the exact solutions to the Schrodinger 
equation for the Gaussian well. For this reason we think that it is interesting to show an 
alternative variational calculation based on trial functions with the expected exponential 
behavior at infinity. 

In Sec. |TT] we develop the variational method for a general single-well potential-energy 
function. In Sec. IIHI we discuss the application of the harmonic and exponential trial func- 
tions to the Gaussian well, compare the approximate and exact energies for the ground 
and first-excited state as well as their shallow-well expansions. In Sec. IIVI we compare the 
variational and exact deep-well expansions that we did not consider in our previous paper.- 
Finally, in Sec. |V]we discuss the main results and draw conclusions. 

II. SIMPLE VARIATIONAL METHOD 

We are interested in the solutions to the Schrodinger equation 

Htp n = e n ip n , n = 0, 1, . . . (1) 

for the dimensionless Hamiltonian operator 

H = f + v(x) = -~+v{x) (2) 

In our previous paper we showed how to convert the Schrodinger equation into a dimen- 
sionless eigenvalue equation and we do not repeat it hereA For concreteness we choose the 
dimensionless potential-energy function to be of the form v(x) = —vof(x), where v o > and 
the function f(x) exhibits a maximum at x = and tends to zero when \x\ — > oo. Since 
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we restrict ourselves to single wells we also assume that f'(x) > if x < and f'(x) < if 
x > although we do not make use of this condition explicitly. We will discuss a particular 
example later on. The reader will find some references about quantum-mechanical wells in 
our earlier paper. 1 

It is our purpose to obtain simple solutions for some states of the quantum well by means 
of the variational method. To this end we choose a variational function ip(a, x) that depends 
on a variational parameter a. If we assume that eo < ei < £2 < • • ■ then the variational 
principle gives us an upper bound to the smallest eigenvalue 

and we realize that the optimum value of a is given by the minimum of W(a). If we write 
W(a) = F(a)-v G(a) where F(a) = (<p\ f \<p) / (<p\ (p) and G(a) = (<p\ f(x) \<p) / (<p\ <p) then 
W'(a) = becomes 

F\a) - v G'(a) = (4) 

Commonly, solving this equation for a may not be possible or may lead to cumbersome 
expressions. In that case we can resort to a simple parametric expression for the approximate 
energy in the form 

F'(a) 



G'(a) 



w{a) = F(a)G'(a) + F'(a)G(a) 

If f(—x) = f(x) then the states have definite parity ip n {—x) = (—l) n tp n (x) and the 
variational principle applies to e or e\ if the variational function is even or odd, respectively. 
Therefore, the variational equations (jSJ) apply to both the ground and first excited states. 

III. THE GAUSSIAN WELL 

In order to illustrate the performance of the variational method on quantum-mechanical 
wells we choose a particularly simple example that allows us to calculate the functions F(a) 
and G(a) analytically. As in our earlier paper we select the Gaussian well given by 

/(*) = e~* 2 (6) 
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and begin present discussion with the calculation based on Harmonic-oscillator-like varia- 
tional functions, i The simplest variational function for the ground state is 



(p oh (a,x) = e ax (7) 



and we easily obtain 



y/2y/a{2a + if' 2 
vo - 2 

/ \ a (4a + 1) 
W 0h (a) = K — '- (8) 

We showed that Woh — > as v o — > according to^ 

W oh = -vl + - 48v% + 832v 8 - 17408w 10 + . . . . (9) 

while the exact expansion is^ 

nvi r- 3 7r(27T + 3v^3 + 3)^ ^^(2^ + 3^)^ 
e — ^- + VZ7ru 1 .... [10) 

Note that the cubic term is missing in the approximate expansion (Q. 
We have also shown that 

<p lh (a,x) = xe~ ax2 (11) 
is a suitable variational function for the first excited state, and in this case we have 

v / 2v / 2aTI (4a 2 + 4a + 1) 

Vo = w* 

W lh {a) = - Q(8a ; + 2a ; 1) (12) 
v ; 2(2a+l) y 1 

It is well known that €\{vq) — > _ as Vq approaches a critical well strength i>o,i ~ 1.342 
from aboveii 2 -. We can estimate the critical strength by means of equations (II 2p in a quite 
simple way. We first obtain a positive root of Wih(a c ) = and then t>o(a c ); the result is 
a c = 1/4 and t>o(l/4) = 9\/3/8 ~ 1.95. Since the variational energy is an upper bound to 
the exact one for all v$ then it is not surprising that v (a c ) > vq \. The error for the critical 
strength obtained with this variational function is rather large: 45%. 

We may try to improve those results by means of more convenient trial functions. Since 
v (x — > ±oo) = we know that the eigenfunctions behave as ip ~ e~ a \ x \ for sufficiently large 



\x\, where a = y/—2e. In order to avoid the function \x\ in our calculations we work on the 
half positive line x > and take into account the boundary conditions at origin: "0(0) ^ 0, 
■0'(O) = for the even states and if)(0) = 0, i/j'(Q) ^ for the odd ones. Thus, the simplest 
trial exponential function for the ground state appears to be 

V9 0e (a, x) = (1 + ax)e' ax , x > (13) 

The calculation of the integrals is straightforward and we obtain 



V ° 0Fe fl2 (4a 6 + 4a 4 - 5a 2 + 2) [1 - erf (a)] - 2a (2a 4 + a 2 - 2) 



a 2 



|v^Fe a2 (4a 6 + a 2 - 2) [erf (a) - 1] + 2a (2a 4 - a 2 + 2) j 



W ° e ^ 10 {0Fe a2 (4a 6 + 4a 4 - 5a 2 + 2) [erf (a) - 1] + 2a (2a 4 + a 2 - 2)} ^ 
where erf (z) is the error function. 

In order to obtain the v -series for W 0e we first expand v (a) into an a-series: 

a 3a 3 4a 4 9a 5 . . 

V0= 2A + V^~^ + ^ + --- (15) 
that we invert to obtain the ^o-series for a 

a = 20^0 - I2^ 2 v 3 + ^—v 4 + 144vr 5 /\ 5 + ... (16) 
Then, we expand Wo e in a Taylor series about a = 

H/ °« = -Io-5 + v5 + '" (17) 

and substitute the series ( 1T61) to obtain the shallow-well expansion 

Woe = _M + Sj^i _ 64^1 + 
5 5 15 

The calculation is straightforward but extremely tedious if carried out by hand. For this 
reason it is a good exercise for showing the students the advantage and power of computer 
algebra. Such software even offer a command for obtaining the inverted series ( TLB"]) in one 
step. For sufficiently small t>o this expansion is slightly better than the one derived in our 
earlier paper, namely Eq. (Q, as follows from comparing the leading terms of the variational 
and exact t> -series: W Q h ~ — Vq, W 0e ~ — 1.26-Uq and e ~ — 1.57vq. However, the cubic term 



that appears in the exact expansion (fTUI) is also missing in the variational treatment based 
on the exponential function f lT3|) . From those results we conclude that e < Wo e < Woh for 
sufficiently small Vq. In other words, the trial function with the correct asymptotic behavior 
at infinity yields a more accurate variational energy for sufficiently shallow wells. 

Fig [TJ shows W h, Wo e and e (obtained by numerical integration 2 ) for some values of v . 
The variational curves are almost indistinguishable in the scale of the Figure. The numerical 
results show that they cross at t>o = v c , where 2.4022 < v c < 2.4023, so that Wo e < Woh if 
v < v c and Wo e > Woh if t>o > v c . In other words: Wo e is more accurate for shallow wells as 
argued above and Woh is more accurate for deep ones. The former inequality is consistent 
with the previous comparison of the leading terms of the fo-expansions. In the next section 
we discuss the deep- well limit. 

For the first excited state we propose the simple trial function 

<p le (a,x) = xe~ ax , x > (19) 

and obtain 



' ° " a [v/?e a2 (4a 4 + 12a 2 + 3) [1 - erf (a)] - 2a (2a 2 + 5)] 



a 2 



jv^Fe" 2 (4a 4 + 8a 2 + 1) [erf (a) - 1] + 2a (2a 2 + 3)} 



Wle< "^ 2{ v /?e a2 (4a 4 +12a 2 + 3)[erf(a)-l]+2a(2a 2 + 5)} (20) 

The variational energy vanishes at a c ~ 0.550 that leads to the approximate critical param- 
eter v (a c ) ~ 1.56 with an error of 17%, quite smaller than that for the harmonic function 
f[TT|) . Once again we appreciate that the exponential function leads to more accurate results 
close to threshold (shallow well). 

Fig. [2] shows Wih, W\ e and t\ (obtained by numerical integration 2 -) for a range of t>o 
values. In this case the variational curves are not so close each other and the crossing takes 
place at 3.5154 < v c < 3.51541. We realize that the exponential function is preferable for 
the description of shallow wells as suggested by the errors in the variational critical well 
strengths. On the other hand, the harmonic function leads to more accurate results for deep 
wells 



6 



IV. THE DEEP- WELL CASE 



The deep- well expansion for a general well can be easily derived by means of perturbation 
theory- It follows from those results that the first terms of the expansion for the Gaussian 
well are 

e n = -v + (n + ^ V^~o - ^ (l + 2n + 2n 2 ) + (y~ 1/2 ^j (21) 

where n = 0, 1, ... is the harmonic-oscillator quantum number. 

In order to obtain the deep-well expansion for Woh first note that the expression for vq 
in equation ([8]) shows that a — > oo as vq — > oo and that the leading term is a w ^p^. 
Therefore, we substitute a = ^^p- + a x + into the expression for v and set the unknown 
coefficients a\ and in order to remove the leading powers of Vq. In this way we obtain 



+ O v~ 6/A (22) 
2 8 1280^ V J K J 

Finally, we substitute this equation into the expression for Woh(a), expand and keep the 

leading terms; the result is 

/ 2^ 3 „ / _i/ 2 



W 0h = -v + ^-- + O [v l/z ) (23) 



Note that this expression agrees with (1211) for n = 0. 

Proceeding exactly in the same way with equation ( !T2|) we obtain 



+ 0^1 (24) 



2 8 128^/vo 
and 

»r u = -„„ + ^ - 15 + o („-v^ (25) 

that agrees with (I2T]) when n = 1. It is clear that the harmonic variational functions are 
suitable for the description of deep wells which is in agreement with the results in figures [1] 
and El 

V. CONCLUSIONS 

We have calculated the first two energy levels of the Gaussian well by means of two types 
of variational functions. The harmonic trial functions lead to simpler expressions and their 
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results are reasonable for all values of the well depth. The approximate energies exhibit a 
slightly accurate expansion for shallow wells and a more satisfactory one for deep wells. In 
the latter case they provide the first three dominant terms of perturbation theory. 

On the other hand, the expressions for the energy obtained from exponential functions 
are rather complicated and their results are less accurate for deep wells. However, they are 
useful if one is interested in the behavior of the system near threshold (shallow wells). 

Present results suggest that the appropriate behavior of the wavefunction at infinity does 
not always guarantee the greatest accuracy. In the case of a deep well the wavefunction is 
strongly localized and compressed about x = and the description of the neighborhood of 
the origin is more important than what happens at larger \x\. 

A most intriguing result is that both variational approaches fail to yield the cubic term 
in the shallow-well expansion. We have carried out a similar calculation for the particle in 
the square box with finite walls and found that the expansion of the variational energy given 
by the trial function ( fl3l) does not exhibit the cubic term that appears in the exact series. 
In this case one can obtain the exact fo-series for eo from the transcendental equation that 
determines the energy levels of that model. We do not have a satisfactory explanation for 
this failure of the variational method. 
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FIG. 1: (Color online) Variational energies Wo e (solid line, blue) and Wo/i (dashed line, red) and 
numerical ones (circles) 




FIG. 2: (Color online) Variational energies W\ e (solid line, blue) and W\h (dashed line, red) and 
numerical ones (circles) 



